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We study a nonadiabatic effect on the conductances in the quantum Hall effect of two- 
dimensional electrons with a periodic potential. We found that the Hall and longitudinal 
conductances oscillate in time with very large frequencies due to quantum fluctuation. 
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Two years after the discovery of the quantum Hall effect, 1 ' 2 Thouless, Kohmoto, Nightm- 
an ' 

gale and den Nijs (TKNN) 3,4 theoretically studied the quantum Hall effect in two-dimensional 

electrons with a periodic potential. They expressed the Hall conductance as an integer multi- 
plied by e 2 /h. The integer is a topological invariant called the Chern number. Such a system 
has been experimentally realized as a superlattice structure in a semi-conductor heterojunc- 

c--: .. ' 

Q tion. 
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We can treat the electric field as a time-dependent vector potential. Then adiabatic ap- 



proximation gives the same expression of the Hall conductance as in the TKNN theory. 8 ' 9 
In this Letter, we consider the conductances of the same system taking a nonadiabatic effect 
into account. We determine the quantum fluctuation of the Hall conductance a xy and the 
longitudinal conductance a yy ; i.e., a xy and a yy oscillate in time with very large frequencies. 

We consider the nonadiabatic effect on the conductances with the help of the Greenwood 
O ■ 

O ■ linear-response theory 10 and study a correction to the adiabatic approximation. Wagner stud- 

ied the failure of the Kubo formula due to a nonadiabatic effect on the dc current in a finite 
system. 11,12 Grimaldi et al. explained the high-T r superconductivity of fullerene compounds 
»■ by taHn g nonadiabatic effects into account* 

This Letter is organized as follows. We first calculate the conductances beyond the adi- 
abatic approximation and obtain correction terms to a xy and a yy of the TKNN theory. We 
then evaluate the oscillations in these correction terms numerically. 

We consider noninteracting electrons in a periodic potential in the x-y plane. A magnetic 
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field is applied in the z direction, and an electric field is applied in the y direction. We treat 
the electric field as a vector potential. Using the Landau gauge, we express the Hamiltonian 
of the system as 

H(t) = ^-(p + eA(t)) 2 + U(x), (1) 

where 

A(t) = (Bx-E y t)e y , 

U(x) = U x cos (—\ + U y cos C^f) • (2) 

We consider the flux per unit cell of U(x) to be the rational number p/q in the unit of 
the flux quantum, which produces p subbands. We label each subband by m (1 < m < p) 
hereafter. Then we obtain the generalized crystal momentum hk defined in the magnetic 
Brillouin zone. 14 < k x < 2ir/qa and < k y < 27r/b. 

We define T~t k (t) as 

H k (t) = e- ikx H(t)e ik - x . (3) 

The instantaneous Hamiltonian TC k (t) is diagonalized as 

Since H k (t) = H kxyky _ eEyt/h {0), we obtain for m / n 

(u mk (x,t)\u nk (x,t)) - 

eE V / du mk( X '°) 



h \ dk y 



ik (x,0))+O(E 2 y ), (5) 



where the dot indicates the time derivative. We note that, by fixing the arbitrary phase in 
each time, we can always achieve (u k (x,t)\u k (x,t)') = 0. 

The time evolution of the density operator is given by the von Neumann equation 

± p (t) = ±[H(t),p(t)}. (6) 

Following the Greenwood linear-response theory, 10 we expand the density operator with re- 
spect to the electric field and take the zeroth- and first-order terms into account 

p(t)=p<® + E vP W(t). (7) 

Thus, the matrix elements 

Pmnk = (^nfc(^)l ^^ P (t)e [kX \u nk (x,t)) (8) 

are approximated using 

- o (0) + E o (1) 

^mnk ^mnk y^mnk 

= JmOmn + E yP mnk ^ V*J 
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where f m (= f(t mk )) is the Fermi distribution. Hereafter, we let e mk denote e mk (0) and 

\ u mk) denote \ u mk( x >°))- 

The time evolution of the off-diagonal elements (m 7^ n) is given by 



d Jl) _ 1 U e U (l) , e / f f J 9u mk 

1P—U - T^( e m k ~ € nk)P, mr ,k + fcUm - Jn) 



dt P mnfc i/j* 



*fc^ft V 



<%. 



:</ 



Therefore, we obtain the matrix elements of the density operator as 



mnt 



du„ 



™*iK 



-ie ■ 



k 



dk„ 



l nk 



Im In 



C mk € nk L 



J_ — Q mr\i nfc 



P (1 \ 

mrnK 



const. 
Let us calculate the current J a (a = x, y) using 

J a = Re[(J a (E y )) - (J a (0))] 



(Ja(E y )) 



Tr p(t) 



dH{t) 



Tr lp k {t) 



8A a (t) 
&H k (t) 



r, -' \ rK "' dk 

By plugging eqs. (11) and (12) into eq. (14), we obtain 



e 2 f d 2 fc v^ 

H Jmbz (2vr) 2 ^ 



<9n 



fc 



<%,, 



ife 



ifc 



<9u 



<% n 



Jfe ) • ( 10 ) 

(m^n), (11) 
(12) 

(13) 
(14) 



(15) 

where MBZ denotes the magnetic Brillouin zone. After some tedious but straightforward 
calculation, we obtain the conductances a xy and a yy in the forms 



a 



•' X 



■<■» 



E ! 

e 



o 



yy 



2irh 



Jy 

E,, 



e 2 f d 2 k v^ 

^Ch + -r TTTTo zZ (fm - fn)r mn (k) sin [(e u - e nk )t/h - mn {k)] , 

n Jmbz (2tt) ^ 



MBZ 



d 2 fc 

(2vr) 2 



E 



du„ 



dk,. 



u 



k 



(fm ~ fn) sin {e mk - e nk )t/h, 



where r mn (k) and 9 mn (k) are the real numbers that satisfy 



r (U\p W rnn(k) 



du 



k 



dk,. 



L nk 



II 



nk 



du 



k 



dk x 



Here, N<jh is the Chern number with a finite-temperature correction 



^Ch = / Jm \ 



d 2 fc 



-Im 



MBZ 7r 



dk,, 



k 



du„ 



k 



dk x 



(16) 



(17) 



(18) 
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We note that if we ignore the time dependence of p . in eq. (10), we obtain a xy = 
(e 2 /2nh)Nch and &yy = 0, which are consistent with the TKNN theory. 

The above-mentioned equations imply that we obtain sinusoidally oscillating terms in 
addition to the Chern-number term in the Hall conductance a xy due to quantum fluctuation. 
The longitudinal conductance a yy also oscillates with zero mean. The oscillation period is 
determined by the energy-gap size. 

Equation (16) is also a gauge invariant. We observe this by treating the electric field as 
the scalar potential Aq{€) = — eE y y exp[— iut] instead of the vector potential. The first-order 
time-dependent perturbation theory gives eqs. (15) and (16) in the limit of 10 — ► 0. 

To demonstrate the oscillation numerically, we assume that the state lies in the lowest 

Landau level. Moreover, we assume that the flux $ (= abB) is given by the unit flux <J>o (= 

2nK/e) multiplied by the rational number p/q as 

g _ abeB _ p 

where p and q are coprime. Due to the periodic potential U(x,y), each Landau level splits 
into p subbands with a g-fold degeneracy in each subband. 

In the weak-potential limit |f7(aj)| <C heB/m c , the perturbation theory gives the instan- 
taneous eigenfunction |« Js)) as 3 

2" 



1 



,k 



x) 



V 



'"%. 



d m (k) 



DO 

'.=— oo 



exp 



\k x I x — £qa 



exp 

uqa 
p 



£ )(-§"'--?)' 



-2my^- 



(20) 



where M is the normalization factor. Here, the coefficients {d 1 ^} and first-order energy shift 
e , are determined by the following secular equation. 

(3d^\k) + a u d v m (k) + (3*d v +\k) = e™. 



d» m {k). 



(21) 



Here, 



n,. 







Uxe -nqb/2pa cog (_ qbky / p + 27TUq/p), 
V -Tvqa/2pb —iqakx /p 



M 



Note that the unperturbed energy e , is independent of m and k; i.e., e r 
After straightforward calculation, we obtain 



'ifl' 



ik e nk 



M) 



(22) 



du mk 


du mk 


dk x 


dk y 


( U nk 


du mk 

dk x 


( U nk 


du u 

rnK 

dky 



V 

£1 


(MiXdd» m 

V dk x J dky 


£«)* 


dk x ' 


v=\ 






V 


r) * 


dd» m 



v=l 



dL, 



ik nk ' 

(23) 
(24) 
(25) 
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Fig. 1. The energy levels around the center of the p (= 65) subbands are shown as functions of (a) 
qak x and (b) bk y . The flux ratio p/q — 65/2 corresponds to B = 10. 4T. We use U x = U y — O.lmeV, 
a = lOOnm, and b = 130nm. The Fermi energy (the dashed line) lies between the 33rd and 34th 
levels. 

Thus, we can calculate the conductances a xy and a yy in eq. (16) using the eigenvalues and 
eigenvectors in eq. (21). 



As typical values, we use a = lOOnm, b = 130nm, and U x 



U„ 



O.lmeV. We assume 



p/q = 65/2, implying B = 10. 4T. With this flux ratio, the lowest Landau level splits into 
65 subbands. We note that the p subbands in the lowest Landau level range from 0.240meV 
to 0.948meV. The center of the next Landau level is located at 1.78meV. In the numerical 
calculation, we obtain the derivatives of d^ik) by the finite-difference method. We consider 
the slices Ak x and Ak y to be 2ir/100qa and 27r/1006, respectively. The Chern number term 
is obtained by Fukui, Hatsugai and Suzuki's method. 15 

First, we make the Fermi energy eF lie between the 33rd and 34th levels, as shown in 
Fig. 1. The fluctuations of a xy and a yy in Fig. 2 consist of the oscillations corresponding to 
the energy gaps between levels higher than ep and levels lower than ep. Figure 2 shows that 
the time resolution required to detect the quantum fluctuation is about lOps. Secondly, we 
make €f lie between the second and third levels, as shown in Fig. 3. Then we obtain almost 
periodic oscillations of a xy and a yy (Fig. 4) in contrast to the seemingly random oscillation 
in the former case. This is because the energy gaps are almost constant in the latter case. 

To summarize, we studied a nonadiabatic effect on the conductances in the TKNN theory. 
We found that both a xy and a yy oscillate due to quantum fluctuation. The oscillation period is 
of the order of lOOps to Ins. It is a challenging but interesting issue to detect these oscillations 
experimentally. 

Although we studied the time dependence of currents at a constant electric field, by similar 
calculation, we observe that constant currents in reverse produce a time-dependent electric 
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Fig. 2. Conductances corresponding to the situation in Fig. 1. In the upper panel, the solid line shows 
a xy as a function of time, and the dashed line shows the Chern number term ach = {g /2irH)Nch- 
In this case, A^ch = 1. The lower panel shows a yy as a function of time. 
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Fig. 3. The energy levels in the lowest part of the p (= 65) subbands are shown as functions of (a) 
qak x and (b) bk y . The flux ratio p/q — 65/2 corresponds to B = 10. 4T. We use U x = U y = O.lmeV, 
a = lOOnm, and b = 130nm. The Fermi energy (the dashed line) lies between the second and third 
levels. In each panel, the lowest solid line below the dashed line expresses the degenerate first and 
second levels. 



field. In this case, the oscillation period is also determined by the energy-gap size. 
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Fig. 4. Conductances corresponding to the situation in Fig. 3. In the upper panel, the solid line shows 
a xy as a function of time, and the dashed line shows the Chern number term ach = {g /2irH)Nch- 
In this case, A^ch = 0. The lower panel shows a yy as a function of time. 
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